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Abstract. The three loop anomalous dimension for the gauge invariant, renormalizable, non- 
local mass operator for a gluon is computed in the MS scheme. In addition the anomalous 
dimensions of the associated localizing ghost fields are also deduced at the same order and it 
is shown that the three loop QCD /3-function correctly emerges from the gluon localizing ghost 
vertex renormalization. 



In non-abelian gauge theories the vector bosons responsible for carrying the quanta of force 
are regarded as massless particles unless there is a spontaneous symmetry breaking. Expressed 
another way there is no gauge invariant local mass operator for gluons in quantum chromody- 
namics (QCD). Whilst it is possible to have BRST invariant masses, such as that which occurs 
in the Curci- Ferrari model, [T], the inclusion of such local mass operators all suffer from the 
disadvantage of leading to non-unitary theories, [21 [3]. Hence they have no predictive power in 
relation to S'-matrix elements. By contrast, there has been an explosion of interest in recent 
years in studying the infrared dynamics of Yang-Mills theories in the infrared limit in the Landau 
gauge using lattice techniques, Dyson-Schwinger equation methods and other more formal ap- 
proaches. One of the main quantities which is analysed is the gluon propagator and it is widely 
acknowledged that it does not satisfy the usual perturbative form of a massless propagator of 
an unconfined field. Instead it is generally fair to say that the gluon propagator, as measured on 
the lattice and other methods, has a behaviour which is not inconsistent with the gluon having 
an effective mass of some sort. Whether this effective mass is due to screening, dynamically 
generated, derived from say Gribov issues, due to vortex condensation or another mechanism 
has not yet been definitively answered. However, if it is to be explained theoretically then one 
is forced into studying extensions of the Yang- Mills or QCD Lagrangians which have a concrete 
gluon mass term of some sort or one where a mass operator condenses. Clearly to do this in 
a gauge invariant way would appear impossible as the obvious mass operator, ^(^4^)^, breaks 
gauge symmetry despite being renormalizable, [T] , where is the gluon field. However, if one 
sacrifices the restriction to local operators then it is possible to have several gauge invariant 
gluon mass terms. In essence there are two types. 

The first, originally introduced in ^ in three dimensions, has been examined in four dimen- 
sions in O E] where it was shown to be renormalizable. Indeed its two loop MS anomalous 
dimension was computed in ^ and shown to be independent of the gauge fixing parameter of a 
linear covariant gauge. The key to demonstrating renormalizability and allowing one to calcu- 
late in a systematic way was the fact that the Lagrangian involving the operator itself could be 
written in terms of local fields additional to the usual gluon, quark and Faddeev-Popov ghost 
fields. These extra (infrared) fields do not affect the usual ultraviolet properties of the original 
non-abelian gauge theory, [5l[6]. Therefore, for example, the /3-function of [7] is unchanged. The 

other type of mass operator is in effect the Stueckelberg term but written as ™? / d'^x [A'^^)^ 
where U is an element of the gauge group. It has been studied in the massive gauge invariant 
model in [HI |9] and is central to a vortex interpretation of confinement. Although also being 
non-local it suffers from the calculational drawback of being non-renormalizable. Though its 
one loop anomalous dimension was computed in [10] for arbitrary linear covariant gauge and 
shown to be independent of the gauge parameter. Part of that calculation rested on the fact 

that the massive gauge invariant operator "^uf J d^x (^^'^)-^ has the renormalizable non-local 
mass operator of [Sj [6] as its first term in a gluon leg expansion of the operator in terms of gauge 
invariant operators, [11]. Therefore the localization of the previous non-local operator into the 
original Yang-Mills fields plus localizing ghost fields provided a useful calculational shortcut. 

From another point of view the non-local operator "^uf J d'^x (^4^^)^ can be viewed as a method 
of gauge fixing QCD in a more concrete fashion as noted in |lllll2t[T3lll4j . This is because that 
gauge fixing operator is gauge invariant and thus avoids the Gribov problem, |15) , which plagues 
the more widely used Landau gauge in the present intense activity into the infrared structure. 

Prom a theoretical point of view one would ultimately like to have a Lagrangian based 
method of studying effective gluon mass which emerges in the current picture and which is 
renormalizable. Moreover, as performing calculations is essential to understanding such low 
energy problems, we focus here on providing the anomalous dimensions of the non-local mass 
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operator of [5,, i6j to three loops in the MS scheme. This is far from being a trivial exercise which 
is due in part to the presence of the additional fields but also because of the generation of a set 
of quartic interactions. As was shown in [5l E] these are essential to preserving multiplicative 
renormalizability. Therefore, we also report on the renormalization of the fields themselves at 
three loops. Indeed as an example of where such three loop results are necessary we note that in 
[161 [IZl the problem of the dynamical generation of a gluon mass was studied in the Landau gauge 
based on the local operator 1(^4^)^. Briefly, the two loop effective potential for this operator was 
computed for Nf massless quarks using the local composite operator (LCO) formalism, \18\ I19j. 
Knowledge of this potential allows one to show that the energetically favoured vacuum is one 
where the operator condenses and therefore dynamically generates a gluon mass. One peculiar 
feature of the LCO formalism, however, is that to have the full two loop potential one needs the 
operator's anomalous dimension at three loops, [161 [171 [IHl [H] . Whilst the results successfully 
demonstrated operator condensation which was stable to loop corrections, [TH [T7|, it suffers 
from one obvious drawback and that is that the calculation was restricted to a specific gauge. 
It would be more appropriate to study the extended operator considered here since it is gauge 
invariant. Indeed this is one of our motivations for this article. However, as will be evident from 
what we present, we believe the determination of this three loop anomalous dimension for the 
non-local operator is sufficiently interesting in its own right to present it separate from an LCO 
computation. 

We begin by recalling the full form of the Lagrangian of [5l [6]. It is 



where a is the linear covariant gauge fixing parameter, is the Faddeev-Popov ghost, ijj'^^ is 
the (massless) quark, i?^^, and H^^^ are the localizing ghosts where the latter are anticommuting 
and m is the gluon mass. For completeness the index ranges are 1 < a < iV^, 1 < / < Np and 
1 < i < Nf where Np and Na are the respective dimensions of the fundamental and adjoint 
representations and Nf is the number of quarks. The covariant derivative, involving the coupling 
constant g, is denoted by D'^ and G^^, is the field strength. The quantities A"'"^'^ are the quartic 
couplings necessary for multiplicative renormalizability and satisfy the symmetry properties 

y^abcd ybacd yabdc y^cdab 

They are not to be confused with the specific rank 4 invariant tensors, such as the totally 
symmetric tensor d'p'^'^ of [2D|, which can be built from the structure functions, /°*'^, or the 
colour group generators, T"". In addition, since the Lagrangian is colour symmetric the quartic 
couplings satisfy a Jacobi style identity, [3 [6] , which is 

j;apq ^pbcd _|_ ^bpqyapcd _j_ jcpq yabpd _|_ jcpqyabcp q 

In ([T|) we have ignored the masses of the {-B^^,, B'^^^jH'^^^, H^^,} sector since they will play no role 
in the present calculation. Finally, we note that ([T|) is the localized version of the Lagrangian 
with the explicit non-local mass operator, 

1 1 - 777^ / 1 \ 

L = - -Gl,G-^- - i^id^Alf - ed^DpC- + i^^'W' - -^Gl, (^-^j G^^^ (4) 
where is the square of the covariant derivative. 
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Clearly with the additional ghost fields and coupling one has to ensure that the gluon, ghost 
and quark anomalous dimensions as well as the usual /3- function remain independent of A"'"^'^. 
This has been verified at two loops in [5j and [6j. Therefore, here we will compute the former 
anomalous dimensions to three loops as well as those for i?^^, and H^^,. The latter will be 
A^^'^'^-dependent. To deduce the anomalous dimension of the mass m or equivalently the mass 
operator anomalous dimension we will renormalize the dimension three gauge invariant operator 
O where 

O = ^ (5;^^ - ^«'^'^) G"^^ (5) 

by inserting it into a gluon B'^^ two point function. The advantage of this approach is that one 
can split the free and interaction Lagrangian in such a way that the operator is in the latter 
and not the former. If it were included in the free part then we would have the huge (and 
unnecessary) computational task of calculating with massive propagators which would require 
the inclusion of the masses of {B^^^, B^,^, H^^, H^j^}. (The explicit mass terms are given in 
[5l[6].) This would be an intractable proposition. Instead treating the operator as an insertion 
means that all fields remain massless and one also avoids the mixing of masses which occurs 
in the full quadratic sector of such a Lagrangian split, aside from the additional complications 
from the 5°^, and H^^^ masses. More crucially with massless fields one can employ the Mincer 
algorithm, [21], which has been encoded, [22], in the symbolic manipulation language Form, [23]. 
The Mincer procedure applies to massless three loop 2-point functions, |2T] , and performs the 
computation in dimensional regularization in d = 4 — 2e dimensions where e is the regularizing 
parameter. Such a high loop order calculation can clearly only be performed via automatic 
Feynman diagram techniques. In such an approach the extraction of the operator anomalous 
dimension is relegated to the evaluation of the divergent part of a 2-point function derived from 
the parent 3-point one, {A'^^OB^^), where the external momentum of the B^^ field is nullified. 
Such a process for this Green's function is infrared safe since no infrared divergent factors such 
as l/(/c^)^ arise in a Feynman integral where k is an internal momentum. 



Green's function 


One loop 


Two loop 


Three loop 


Total 


Aa Ab 


5 


52 


1279 


1336 




1 


8 


152 


161 




1 


8 


152 


161 


B" B^ 

flu crp 
jra Trb 

A^pBlB^p^ 
KOBl 


1 


20 


464 


485 


1 

7 
5 


20 
166 
131 


464 
5827 
6917 


485 
6000 
7053 


Total 


21 


405 


15255 


15681 



Table 1. Number of Feynman diagrams for each Green's function. 



In the final part of this setup description we note that we have to be careful in ensuring 
the correctness of the final expression. Since the operator insertion is in a Green's function 
involving a localizing ghost, we require a strong check on the -B^^ renormalization constants. 
To ensure this we have also performed the three loop MS renormalization of the A'^B^^B'^^^ 
vertex itself. As in the original QCD Lagrangian, this Green's function will produce the three 
loop MS /5-function of the gauge coupling, [71 [Ml US]. Again as this is a 3-point function we 
nullify the external momentum of the B'^^ field relegating it to a 2-point function whence it can 
be determined by the Mincer algorithm. In Table 1, we have listed the number of Feynman 
diagrams computed for the present article. Those for the gluon and ghost exceed the numbers 
for the corresponding original QCD calculations due to the presence of the localizing fields. The 
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numbers of graphs in Table 1 are deduced from the Qgraf package, [26 j , which is the starting 
point for each of the Green's functions. The Qgraf routine generates the Feynman diagrams 
electronically and these are then converted to Form input notation prior to the application of 
the Mincer algorithm. One additional complication is the non-trivial task of extending the 
Form group theory module to handle the group theory associated with the A"^^*^ couplings 
subject to the symmetry and Jacobi properties of ([2]) and ([3]). In addition we have also used the 
property noted in [6] that 



y^acde yhcde 



^ab ycdpq \^cdpq 

Na 



yacde yhdce 



^ ^ab \cdpq \cpdq 

Na 



(6) 



and the analogous extension to the products of three A'^'^^'^-tensors with two free indices, which 
follow from the fact that there is only one rank two isotropic tensor in a classical Lie group. 
Finally, we note that the propagators of the (massless) localizing ghosts are given in 



We now record our main results at three loops. First, we define the renormalization constants 
for the relevant fields and the operator as 



VZbB"'"' , H^f"" = VZhH^'"' , Co = ZoO (7) 
where the subscript o denotes the bare quantity. Then the respective anomalous dimensions are 



7B(a, A, a) 
where we note 
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70 (a. A) 



a 
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ail 



oa 



(8) 



(9) 



with j3{a) the /3-function of the gauge coupling a = (7^/(167r^) and 7^(0,0) is the anomalous 
dimension of the linear covariant gauge fixing parameter. We use the conventions of [27] for this 
and note that 7a (a, a) = — 7^(0, ct) with the latter defined to be the gluon anomalous dimension. 
In the term involving the /3-function of the quartic couplings, /3"^'^'^(a, A), it is understood that 
the differentiation respects the symmetries of the A^^'^'^-tensor couplings, ([2]). Given this we find 
the MS expression 



7B(a, A, a) 
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(a-3)CA + 
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rf^yabcd ^^acpq ybpdq _j_ yabcd ^^apcq ybqdp 



+ 0{aS\' 



(10) 



where C{n) is the Riemann zeta function and a factor of l/(47r), which derives from the loop 
integral measure, has formally been absorbed into each A"^'^'^ to simplify the presentation. The 
result ([Toll explicitly verifies the equality of the Slavnov-Taylor identity of [5l [6] to three loops. 
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The group Casimirs are defined by TV (t'^T*) = TfS''\ T'^T" = Cf and p^'^f'^d ^ (j^^a 
have introduced the shorthand notation for the contraction of two structure functions 



rabed 
JA 



jeab jeed 



'. We 
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and defined the order symbol, 0(a^;A^), to correspond to the four loop corrections. Moreover, 
with these values the three loop MS QCD /3-function of [3 [2^ [25] correctly emerges as A'^*'^'^ 
and a independent from the A^B^^B'^^ vertex. Given this we find the three loop correction to 
the gluon mass operator is 



1 



3 

abed ^aebd 



a + 
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24' 
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768" 



fl2) 



which is clearly a independent as expected on general grounds but which in fact provides a non- 
trivial check on our computation. It is worth stressing that the emergence of the correct A"^*^*^ 
independent /3-function and the gauge parameter independent 70(0, A) is a non-trivial check on 
the implementation of the symmetry properties of A"'"^'^ in the FORM group theory module. 

One additional calculational detail is worth noting and that is that A"^'^'^ itself undergoes a 
renormalization within the three loop calculations. Its one loop /3-function was given in [6] as 



nabcd 



(a, A) 



i(d-4)A 



abed 



1 

+ 



yabpq ycqdp _j_ ^apbq ^^edpq _|_ yapcqybpdq _j_ yapdq ybpcq 



apbq nepdq 2 



(13) 



Since we are going one loop beyond [6], it might have been expected that the two loop MS 
correction of (jl3p was needed. However, for the operator renormalization the first place A"^^'^ 
occurs is at two loops. Therefore, one only needs its one loop renormalization. Equally for the 
B'^j^ and H^^^ anomalous dimensions A"*'^'^ first appears at two loops and again only its one loop 
renormalization is necessary to deduce the fully renormalized three loop 2-point function. Here 
this is because the one loop graph involving A"'"^'^ which contributes to either 2-point function 
results in a snail graph which clearly is zero for the massless fields we consider. In other words 
it would only contribute to the renormalization of the B^^, or H^^^ mass renormalization. Such a 
property of the A^^'^'^-structure of the anomalous dimensions in fact prevents us from having to 
extend the A"*'^'^ renormalization to two loops by renormalizing massless 4-point functions which 
have non-safe nullifiable external momenta and hence not accessible to the Mincer algorithm. 
Moreover, it is worth noting that this is the first use of ()13p within a computation and the overall 
consistency of our three loop renormalization is a non-trivial check on its correctness. 

As a final check on our anomalous dimensions, we note that in the original renormalization 
constants we have been careful to check that the triple and double poles in e are correctly 
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predicted from the known one and two loop structure. For the current Lagrangian, ([T]), this has 
an additional feature and that is that one has to take into account two coupling constants, a 
and X"-^'^'^, To aid the interested reader in this respect, we provide the explicit three loop MS 
renormalization constant for O whence ()12p was deduced. It i^ 
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We note that the extraction of all our renormalization constants followed the procedure discussed 
in [28] for automatic Feynman diagram computations where all the Green's functions are first 
determined as a function of the bare quantities, such as the coupling constant. The renormalized 
versions are then introduced by the standard rescaling and the remaining divergences in e are 
absorbed into the appropriate renormalization constant associated with that particular Green's 
function. This procedure is straightforward to implement within the Form routines. 

To conclude we have provided the full three loop MS renormalization of all quantities in 
the localizing Lagrangian, ([1]), for a renormalizable gauge invariant non-local mass operator, 
except the quartic couplings. Whilst this has been a cumbersome task it will actually play 
a crucial role when one considers the corresponding renormalization of the full operator of 

"{c/f / d^x (A^'^)2 beyond the one loop result of [TD]. Although this is clearly non-renormalizable, 
it is not inconceivable that one could localize the operator order by order in perturbation theory 
with a finite set of localizing fields at each order. Indeed the one loop result of [lOj proceeded 



^We note that attached to the version of this article which appears on the arXiv there is a Form file which 
contains the results (fTO)l . (fT2)) and (fTi)) . 



7 



under this assumption. Whilst the additional couplings are absent at that level it would be 
interesting to see the structure of the gauge parameter independent anomalous dimension which 
emerges and to study the role of extra quartic couplings play in any renormalization group 
evolution. Moreover, given the successful extraction of the three loop anomalous dimension, it 
now also opens up the possibility of computing the two loop effective potential of this gauge 
invariant operator to study its condensation properties. 
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